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We study correlators of the global U{1) currents in holographic models which involve 
J\f = 4: SYM coupled to the finite density matter in the probe brane sector. We find the 
spectral density associated with the longitudinal response to be exhausted by the zero sound 
pole and argue that this could be consistent with the behavior of Fermi liquid with vanishing 
Fermi velocity. However the transversal response shows an unusual momentum independent 
behavior. Inclusion of magnetic field leads to a gap in the dispersion relation for the zero 
sound mode propagating in the plane of magnetic field. For small values of the magnetic field 
B the gap in the spectrum scales linearly with B, which is consistent with Kohn's theorem 
for nonrelativistic fermions with pairwise interaction. We do not find signatures of multiple 
Landau levels expected in Landau Fermi liquid theory. We also consider the influence of 
generic higher derivative corrections on the form of the spectral function. 



I. INTRODUCTION AND SUMMARY 



Perhaps the deepest open problem in condensed matter physics is the classification of compress- 
ible quantum liquids. This refers to stable states of zero temperature quantum matter that do not 
break any symmetry and support massless excitations. This question cannot be easily addressed 
within the confines of standard field theory. The issue arises when fermions are considered at 
finite density and the culprit is known as the "fermion sign" problem. Dealing with time-reversal 
symmetric finite density bosonic matter the methods of equilibrium statistical physics give a full 
control and invariably one finds that the ground states break symmetry. Dealing with incompress- 
ible quantum fluids like the fractional quantum Hall states the mass gap is quite instrumental 
to control the theory, revealing the profound non-classical phenomenon of topological order. The 
hardship is with the compressible quantum fluids: the only example which is fully understood is 
the Fermi-liquid. 

The ease of the mathematical description of the Fermi-liquid as the adiabatic continuation 
of the Fermi-gas is in a way deceptive. Compared to classical fluids its low energy spectrum of 
non-charged excitations is amazingly rich. In addition to the zero sound, there is a continuum 
of volume conserving "shape fluctuations" of the Fermi-surface, corresponding with the particle- 
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hole excitations (Lindhard continuum) of the conventional perturbative lore. Although serious 
doubts exist regarding the mathematical consistency and their relevance towards real physics, the 
"fractionalized (spin) liquids" that were constructed in condensed matter physics appear to be 
still controlled by the presence of a Fermi-surface while these are not Landau Fermi-liquids in the 
strict sense. This inspired Sachdev to put forward the interesting conjecture that the Fermi-surface 
might be ubiquitous for all compressible quantum liquids [Ij. 

The gauge-gravity duality or AdS/CFT correspondence provides a unique framework to deal 
with these matters in a controlled way (see [l-6j for recent reviews). Although it addresses field 
theories that are at first sight very remote from the interacting electrons of condensed matter, 
there are reasons to believe that it reveals generic emergence phenomena associated with strongly 
interacting quantum systems. Field theories whose understanding is plagued by the "fermion 
sign" problem appear to be quite tractable in the dual gravitational description. With regard to 
unconventional Fermion physics, perhaps the most important achievement has been the discovery 
of the "AdS2 metal" PH'S], dual to the asymptotically AdS Reissner-Nordstrom black hole. On 
the field theory side this describes a local (purely temporal) quantum critical state that was not 
expected on basis of conventional field theoretic means. Although quite promising regarding the 
intermediate temperature physics (the "strange" normal states) in high Tc superconductors and 
so forth, this AdS2 metal is probably not a stable state, given its zero temperature entropy. Much 
of recent activity has been devoted to the study of the instability of this metal towards bosonic 
symmetry breaking (holographic superconductivity ^Oj, "stripe" instabilities [10]) and towards the 
stable Fermi liquid (llH13j . 

The top-down constructions might become quite instrumental in facilitating the search for 
truly new quantum liquids. An important category are the Dp/Dq brane intersections; the p = 3 
case provides us with a set of especially tractable examples. The dynamics of the low energy 
degrees of freedom of the D3-Dp strings can be studied in the probe approximation where the 
back- reaction to the AdS^ x geometry can be neglected IH]. In this paper we will consider 
D3 and Dp branes intersecting along 2+1 dimensions, where p=5 (p=7) corresponds to the (non)- 
supersymmetric system. As emphasized in [15J the nonsupersymmetric system can be viewed as 
a model of graphene: the brane intersection fermions are like the Dirac fermions moving on the 
2+ ID graphene backbone, (tunable to finite density by gating), interacting strongly through the 
gauge fields living in 3+1 dimensions. We will present a number of results for the longitudinal- 
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and transversal dynamical charge susceptibilities (at finite frequency w and momenta q)^, in the 
absence and presence of a magnetic field, for both the supersymmetric and non-supersymmetric 
D3/Dp systems at finite density. We find very similar results in both the supersymmetric- and 
fermionic set ups, showing that these outcomes at strong 't Hooft coupling are not caused by the 
difference in the Lagrangians. We find suggestive indications for the presence of an entirely new 
form of quantum liquid, but we cannot be entirely conclusive. Our observations cannot entirely 
rule out the existence of a Fermi liquid with vanishing Fermi velocity. 

In fact, the first study of these systems at finite density already produced evidence that some 
odd state is created. In ref. [TB] it was observed that the density-dependent part of the heat 
capacity in the D3/Dp systems with 2 + 1 dimensional intersection behaves like T^. This is in 
contrast to the result for the Fermi-liquids which is set by the Sommerfeld law of the specific 
heat C = ^T, where the Sommerfeld coefficient 7 is proportional to the quasiparticle mass. This 
behavior remains to be understood: for example, it is conceivable that the linear term in the heat 
capacity exists, but is parametrically suppressed in the holographic model. On a side, it is worth 
noting that in the context of pnictide superconductivity a rogue signal has been detected that 
refuses to disappear: this indicates that the electronic specific heat of the metal state ~ |17j . 

As mentioned above, besides the Lindhard continuum an interacting Fermi liquid will carry a 
single propagating mode called zero sound. Unlike the usual sound at finite temperature, transla- 
tional invariance alone is not sufficient for establishing the existence of the zero sound mode. The 
discovery of zero sound associated with the brane intersection matter |l6j is therefore significant. 
The fate of the holographic zero sound was further studied in p^H25] (see also |26l ETj for closely 
related work) . At very low temperature the attenuation (damping) of this zero sound behaves like 
the ( "collisionless" ) Fermi liquid zero sound, in the sense that it increases like the square of its 
momentum. In |24j it was found also that upon increasing temperature the zero sound velocity 
decreases while the attenuation increases, turning into a purely diffusive pole at high temperatures. 
This is different from the crossover from zero sound to ordinary sound as function of temperature 
in a single component Fermi-system like ^He. In the brane intersection systems momentum is 
shared between the superconformal strongly coupled uncharged sector and the material system on 
the intersection, and the latter does not support hydrodynamical sound in isolation. Somehow, 
upon lowering temperature the momentum of the brane intersection matter becomes separately 



^ In this paper we denote the values of frequency and momentum by bold letters. The usual letters, defined below, 
are reserved for dimensionless variables. 



4 



conserved, facilitating the emergence of the zero sound in the low temperature limit. 

Given that zero sound is rather ubiquitous, one would like to obtain more direct information 
regarding the density fluctuations of the quantum liquid. These are expected to be contained 
in the fully dynamical, momentum and energy dependent charge susceptibility/density-density 
propagator associated with the conserved charge on the brane-intersection. One strategy is to look 
for the momentum dependence of the reactive response (real part) at zero frequency: one expects 
a singularity at twice the Fermi momentum, 2qp where the Luttinger's theorem implies that qp is 
set by the bare chemical potential, ~ . A number of papers has been devoted to the search 
of such singular behavior in the framework of AdS/CFT. In [12] the {J^J^) correlator has been 
computed in the holographic setup where the only charged degrees of freedom are four-dimensional 
fermions. The resulting function was completely smooth. In |28H31j the two-point function for 
global currents was computed for various systems and again the tree-level computation in the bulk 
did not show any nonanalytic behavior. Very recently it has been argued that a singularity can be 
observed in the systems where an exact result to all orders in a' is available [32]. 

Searching for the singularity at 2qF is in principle a tricky procedure because these "Friedel 
oscillation" singularities are strongly weakened by the self energy effects in the strongly interacting 
Fermi-liquid. Another way to probe for the signatures of the Fermi liquid is to compute the 
imaginary part of the dynamical density susceptibility in a large kinematical window because this 
spectral function shows directly the density excitations of the system. The result is well known in 
the weakly interacting Fermi liquid, see Fig. 1: besides the zero sound pole one finds the Lindhard 
continuum of particle hole excitations. It is worth noting that as the value of the Landau parameter 
Fq increases, the spectral weight in the density response is increasingly concentrated in the zero 
sound poles, "hiding" the Lindhard continuum. In this regard the transversal density propagator 
is quite informative: since in this channel no collective modes are expected to form, this is the 
place to look for the incoherent Fermi-surface fluctuations. Unfortunately technical issues prevent 
us from accessing the regime of parametrically small Fermi velocity. Our holographic computations 
of the longitudinal- and transversal dynamical charge susceptibilities are limited to a kinematical 
window where w ~ |q|. 

Despite this caveat, the holographic density propagators that we compute reveal very interesting 
information. We find that the longitudinal density propagator is within our numerical resolution 
completely exhausted by the zero sound pole (Fig. 4). Regardless the precise nature of the un- 
derlying state this signals very strong density /density correlations in this liquid. The transversal 
charge propagator shows that sound is not the whole story. The "other stuff', albeit very un- 
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like a Lindhard continuum, signals the presence of a sector of highly collective, deep IR density 
fluctuations: the imaginary part of the transversal propagator behaves like x|*''(q, w) ~ w. This 
response is surprisingly momentum independent and suggests local quantum criticality, which was 
instrumental in the "AdS2 metal" setup. All of this seems to imply that we are indeed dealing 
with some entirely new quantum liquid. 

To probe some of the features of this quantum liquid, we introduce an external magnetic field 
which is a valuable "experimental tool". This induces the gap in the spectrum that is visible in 
the holographic calculations. Dealing with a 2+lD Fermi-liquid one would expect the signatures 
of Landau levels also in the density response. In the strongly interacting system, the longitudinal 
response should reveal the "magneto-roton" , the left over of zero sound in the system with a 
magnetic field which is well known from (fractional) quantum Hall systems [33j ^. According to 
Kohn's theorem [35j, the density spectrum should show a gap equal to the cyclotron frequency at 
zero momentum. Note that this theorem is very generic and only assumes that degrees of freedom, 
charged under the magnetic field, interact pairwise. Our holographic calculation reveals that: i) 
at small values of the magnetic field B the value of the gap'^ scales linearly with B, which is 
consistent with Kohn's theorem for the nonrelativistic fermions and ii) there are no signatures of 
Landau levels associated with incoherent particle- hole excitations (Fig. 2). 

The remainder of this paper is organized as follows. The next section is devoted to the review of 
Landau Fermi liquid theory including the random phase approximation (RPA) for the dynamical 
response. In particular, we review the appearance of the zero sound mode in the RPA calculation 
of the density- density correlator. As the value of the interaction strength increases, the Lindhard 
continuum gets separated from the zero sound pole (Fig. 1) and gradually disappears. In the 
extreme limit of vanishing Fermi velocity, the spectral density is completely exhausted by the zero 
sound mode. We also review the RPA expectations for the 2+1 dimensional fermion system in the 
presence of magnetic field. There we expect Landau levels to contribute to the spectral density 
(Fig. 2). 

In Section III we review the holographic description of the D3/Dp brane systems. The subject 
of our interest is the fermion matter, which is formed (at finite chemical potential for the fermion 
number) in the low energy theory living on intersection of the Nc D3 branes and Nf Dp branes. 
We consider the case of Nc ^ Nf ~ 1 and strong 't Hooft coupling A, where the holographic 

See [35 for related work in the context of holography. 
^ This is also consistent with the observations made in [231 136| where the same D3/D7 system, modified by the 
inclusion of flux through the internal cycles, is considered. 
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description is applicable. 

In Section IV we focus on the zero sound mode and show that it develops a gap in the presence 
of magnetic field. In the case of vanishing magnetic field, B = 0, we observe a zero sound mode 
whose speed is the same as that of the first sound. As long as the value of the magnetic field B 
is small compared to w^,q^ (in appropriate units), the sound mode peak in the spectral function 
is not significantly affected. On the other hand, the presence of the nonvanishing magnetic field 
leads to a gap in the dispersion relation for zero sound. (The effective action proposed by Nickel 
and Son [37] in the presence of the magnetic field gives vanishing sound velocity). In the regime 
of small magnetic field we derive the scaling behavior of the gap in the spectrum Wc as a function 
of magnetic field. The result, Wc ~ i? is consistent with fermions acquiring an effective mass. 

In Section V we investigate the current-current correlator at non-vanishing frequency w and 
momentum q. We observe that in the longitudinal channel, the only nontrivial structure both 
in the real and in imaginary parts of the correlators is provided by the zero sound. There is no 
nontrivial structure in the transverse correlators when B = 0. We discuss our results in Section 
VI. 

In Appendix we consider higher derivative corrections and show that when they are added to 
the DBI the correlators are not significantly modified. 



II. FERMI LIQUID AND THE RANDOM PHASE APPROXIMATION 

In this section we review the application of the random phase approximation (RPA) for the 
computation of the density- density response function ( Jo(w, q) Jo(w, — q)) in Landau Fermi liquid 
theory. We consider the 2-1-1 dimensional theory for both cases of vanishing and non- vanishing 
magnetic field. 

Due to the interaction of quasiparticles, the variation of quasiparticle energy due to small 
perturbation of the distribution function, is given by (see, e.g, [38]) 

fe(q) = J dq7(q,q)'^"(q') (1) 

Because the small changes of quasiparticle density occur in the vicinity of a Fermi surface, one 
considers the function /(q, q') to be dependent on the momenta on the Fermi surface, and therefore 
it boils down to a function of the angle between q and q': 

777 

—m = 2F{e) . (2) 
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where, as usual, the effective mass at the Fermi surface is defined via 

* qF <9e(q) 

m = , VF = q=qj. (3) 

Landau parameters Fi are the coefficients of the expansion of F{6) in Legendre polynomials: 

F{e) = J2i'^l + l)FiPiicose) (4) 

The Fermi liquid has a collective excitation at vanishing temperature called zero sound. In the 
case of i^; = 0, I > 0, the speed of zero sound uq can be determined from 

s s + 1 1 uo 

-log --1 = — , s= — (5) 

2 s — 1 Fq vp 

which, in the limit i^o S> 1 gives s ~ ^/Fq- 

To compute the dynamical collective responses of a Fermi liquid, one evaluates the time depen- 
dent mean field (random phase approximation) obtained by summing up the quasiparticle "bubble" 
diagrams. Assuming for simplicity only the presence of a contact interaction, with effective cou- 
pling constant V ~ Fq, the nth diagram is equal to y"~-'^(xo(q) w))". The susceptibility in the 
RPA is then given by the sum of a geometric progression: 

I . Xo(q,w) 
x(q,w =- — -, (6) 

i-^xo(q,w) 

Express x = X^*"^ + ^X^*\ hence 

{i-Vxt\^,^)? + {xf{^,^)?' 

Then we study density of excitations by plotting x^^\^.-,^)- The result for vanishing magnetic 
field is presented in Fig. [T| where we plot the susceptibility (for qp = 0.2) at strong and weak 
coupling V . In the case of strong coupling there is a finite gap, separating the zero sound collective 
mode, and the band of the particle-hole excitations. For given small frequency w the width of the 
gap is given by (5q ~ ~ Note the non-analytic step behavior at q = 2qi7, originating from 
the free response function Xq\c{,w). In the case of weak coupling the zero sound mode merges 
with the left edge of particle-hole band. 

The location of zero sound pole is determined as a solution to equations Xo^Cq^w) = 0, 
Xo^^(q, w) = 1/V. The real part Xo^'^CqiW) of Lindhard function for 2D Fermi gas is given by 
(see, e.g., [IQ]): 



Xr(q,w) = - 1 + 



sign(i/_)0(|z/_| - l)Jvl - 1 - s\gn{iy+)e{\u+\ - l)Ju^ 



(8) 



FIG. 1: Spectral density x'-^HQi'^) strong coupling (V = 50, left graph) and weak coupling {V — 3, 
right graph) in the random phase approximation, at vanishing magnetic field. Fermi momentum is put to 
qp = 0.2. Note that at strong coupling zero sound is well separated from the particle- hole continuum, while 
at weak coupling zero sound merges with the left edge of the particle-hole continuum. At small frequencies 
particle-hole continuum sharply ends at q = 2qp. 

where i'± = For large = s ^ 1 one may expand 

X^^(q,w) ~ ^-f . (9) 

Therefore, for the speed of zero sound one obtains s = \JV jl^ exactly as it follows at large Fq from 
the equation ([5]). 

Suppose now that besides Fq there is also non- vanishing "mass" Landau parameter F\. In the 
relativistic case, the value of m* is related to the value of the chemical potential [39j, 

"^* = ^(l + y) (10) 
The speed of zero sound then satisfies equation 



For free fermions in a magnetic field B, the Lindhard function is equal to (see, e.g., |40j ) 

Xo(q,w = — 2> — — q) , (12) 
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where 



(13) 



for n' > n. Here we have introduced the cyclotron frequency Wc = B/m* and the magnetic length 
I = The functions L'^~'"' are Laguerre polynomials, and f{£n) is an occupation number for 

the nth Landau level. 

We would like to compute the effect of the magnetic field on the density-density response 
function of the interacting fermions. Let us write the quasiparticle interaction Hamiltonian 

Hint = ^ Vqnqn_q (14) 

q 

in the basis of Landau levels wavefunctions. The corresponding matrix elements of the density 
fluctuation operator nq = '^ii'^t.+q given by 

/ qx(ky + k')£2\ 

(n ky\nci\nky) = exp I -t I F„/„(q)5ky_k;,q., • (15) 

The density fluctuation operator in the basis of Landau level wavefunctions is then given by 

Note that 

{{nky\nq\nky)y = (n'ky|n_q|nky) (17) 



implies (n-q)^ = n_q. Substituting (16) into the interaction Hamiltonian (14), assuming again only 
a contact interaction of plane waves Vq = V ^ Fq, and considering all quasiparticles in the same 
Landau level n, one obtains 

H^nt = V Cnk,4^_q/nk^4, +q,^ exp (-^e\A^y - - q,) - ^) [L^M^i' ■ ^ 

Let us choose the momentum to be in y-direction, then 

Hint = ^ y^yCn'kyC^nky-Ciy^nk.'yC^nk'y+ciy ' (1^) 
'iyM-yM-'y 



where ^q^ = [L^(q^£V2)]^ exp [-^) V. 

We can explicitly demonstrate that the zero sound mode is gapped in the magnetic field, with 
the gap being equal to Wc, in agreement with the Kohn's theorem [35j. For this aim we are to 
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solve equation Xo'\it'^) = 1/^ again. From (12), (13) one may obtain the following expression 
for Xo : 



k=l j 



{j + k)\ 



— (A;wc)-^ ' 



(20) 



where the prime denotes summation in the range max(0, f — k) < j < v, and v is the number of 
occupied Landau levels. Following |40j . we consider this equation for small q and w ~ Wc. Then 
the main contribution in the sum over k comes from the term with k = 1, and we obtain equation: 



const - 



1 



— V ' 

and therefore the zero sound dispersion relation is given by 



(21) 



w 



(22) 



where c ~ Vwc is a constant. Similarly, for any integer M, there is a mode with dispersion relation 



w = J(Mwc)2 + c'q2A^. 



(23) 



We plot RPA computations of two-point function, for Uc = 0.25, restricting to the first two first 
branches, in Fig. [2j 




FIG. 2: Spectral density in the random phase approximation of the 2 + 1 dimensional Fermi liquid in the 
plane of the magnetic field, with cyclotron frequency uJc = 0.25. First two of infinitely many collective 
excitation branches are shown. Each branch starts at {q — 0, uj — Mujc), where M is an integer. 
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III. Dp BRANE IN AdS^ x BACKGROUND 

We study strongly interacting massless fermions at zero temperature and finite density. A good 
field theoretical model of such a system is A/" = 4 SYM theory with gauge group SU{Nc), coupled 
to matter in the fundamental representation. A convenient way to study strongly coupled theories 
is provided by holography where one considers a dual gravitational theory, taking the limit of large 
't Hooft coupling A = gyi^^jNc, and the limit of large Nc- The dual gravitational background is 
created hy Nc ^ 1 D3 branes, and has an AdS^ x geometry. The coupling to fundamental 
matter is realized by considering an embedding of a probe Dp brane in the AdS^ x background 
|14j . We will consider D3/Dp configurations with d = 2 + 1 dimensional intersections. 

Let us now provide a more detailed description of the bulk gravitational theory set-up. Consider 
AdS^ X geometry, with the metric 

ds^ = I r'^{-dt^ + + ^ + dOf J . (24) 

Here L is the radius of and scale of curvature of ^^5*5. We will study the probe Dp brane. 



embedded in the geometry described by ( 24 ) . We represent the metric on as 



d^l = dnl + sin^ 9dnl_^ = d9^ + sin^ 9dnl_^ + cos^ 9dnl_^ , 
where n = p + 1 — d. Then we define coordinates p , / via the relation 

p = r sin 61 , f = rcos9 , r"^ = + f , (25) 

and write 

d9^=^-l^J^dp\ dr^=^-P±lp^dp\ (26) 
which gives the following induced Dp brane world-volume metric 

d4p = (^r\-dt^ + dxidx') + ^ (1 + (dpff) dp^ + ^df^Li) • (27) 



The coordinate f{p) defines an embedding of the Dp brane in the AdS background (24). In the 
case of the trivial embedding f{p) = 0, which is what we are going to deal with in this paper. Dp 
brane crosses the Poincare horizon of the AdS space. In the case of d = 3 p = 7 such a configuration 
becomes stable only for sufficiently large values of chemical potential pch in the dual field theory 
|41j . (See also |42j for the phase structure of the similar model in the presence of the magnetic 
field.) Note that holographically computed correlators do not depend on the dimensionality of the 
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probe brane; in particular our results apply in the case of stable supersymmetric D3/D5 defect 
theory. 

Subsequently we add a gauge field on the world-volume of the probe D7 brane. In general 
we are interested in non- vanishing magnetic field B. So we consider the following components of 
the field strength: 

Fi2 = 5, Fop = -apAo(p) . (28) 
Consequently the DBI action for the Dp brane is given by ^ 

Sdbi -Jjj dP+'x^/-det{G + F) = j dn^-i j d'^xS, (29) 
where we have denoted 

S ~ iVei/-5 j dpp''-^ ^ {L^p^ + B^){1 - {dpAofL-^) . (30) 
Now rescale gauge field on the world-volume as 

^, = ^, (31) 

which yields the DBI action in the form, 

S ~ N,LP-^ J dpp'^-y(p4 + 52)(l-(9pAo)2), (32) 

where B = B/L'^. 

In the case of a non-vanishing magnetic field there is also a Chern-Simons term in the total 
action for the Dp brane. It can be shown that this term vanishes in the case of / = embedding. 

The boundary value of Aq is equal to the chemical potential of the dual field theory: Ao(p = 
oo) = P-ch- Due to f{p = 0) = and the initial condition AQ{r = 0) = (imposed to ensure 
that chemical potential vanishes when the charge density is zero) we obtain Ao{p = 0) = 0, and 
therefore the chemical potential may be expressed as 



oo 



IJ-ch = / dpdpAo . (33) 

^0 

Introducing a constant of integration d, the solution of the equation of motion for dpAo field 
strength becomes, 

dpAo = . (34) 



We adopt the convention 2ira' = 1. For our purposes we are ignoring the total numerical coefficient, which leaves 

9s ~ A ~ pt 



US with an overall normalization of the action proportional to — - - . 
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Using this expression and eq. (33), we obtain the value of the chemical potential 

4r(5/4)2 J2 



^ch = / dp dpAo 

10 



(35) 



IV. HOLOGRAPHIC ZERO SOUND 

In this and the next sections we study D3/Dp system with d = 2 + 1 dimensional intersection, 
described by trivial /(p) = embedding of the probe Dp brane in the AdS^ x background. We 
consider the gauge field on the Dp brane world-volume, solve its classical equations of motion and 
use AdS/CFT to find the two-point functions of the U{1) current in the dual field theory. In this 
section we show the existence of holographic zero sound in the D3/Dp configuration, to observe 
that it develops a gap as the magnetic field is turned on. In the next section we will study the 
current-current correlation function numerically. 

A. Zero sound in the D3/Dp system with d = 2 + 1 dimensional intersection 



Equation (34) is the expression for the background field strength dpAo. Let us turn on small 



fluctuations oq, ai, 02, dependent on coordinates x^, x^, p. In addition let us fix the gauge dp = 0. 
The longitudinal response is described holographically by the ao and 0,2 components of the gauge 
fleld, and the transverse response is described by the ai component. The DBI action, expanded 
up to the second order in fluctuations, then takes the form ^ 

(dpaof 



S 



dp 



1 - {dpAoy V 1 - (dpAof 



^ P^{dpa2f - {doa2 - ^200)^ ' _^ 



p4 + ^2 



+ 



+ 



p^ + B^ V / + 
2BdpAo 

(p4 + S2)(l-(9pAo)2) 



1 - {dpAo^ 
{d2aidpao - 9oai5pa2 + (5oa2 - d2do)dpdi) 



(36) 



Note that the last line in ( 36 ) describes a coupling of the transverse and longitudinal gauge potential 



components. Bellow we will consider Fourier transform of the gauge field 



(2vr)- 



-e-*^"°+^'i"^a^(p,w,q) 



(37) 



^ We thank J. Shock for comments on this action. 
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Now we substitute eq. (34) into the action (36), define 6^ = S^/d^, and introduce a new variable 



z = j^, so that z = is a boundary and z = oo is a Poincare horizon of AdS^. In addition, we 
make the quantities w, q dimensionless, by measuring these in units of d: w = ojd, q = qd. We 
also denote for shortness of notation 



C = 1 + {1 + b')z 



S 



(38) 



26z''((?2aiaQ - 9oai«2 + (c^o«2 - 52ao)a'i)) , (39) 



Then the action ( 36 ) becomes written as 

dz 



1 + 62^4 



where we have omitted bars for simplicity of notation, and prime denotes differentiation w.r.t. z. 
In momentum representation 



S 



dz 



1 + 62^^ 



(-C'/2^[,(w, g)a'o(-w, -q) + Q^'^a'^{uj, q)a2{-co, -q) + C^/^E{u;, q)E{-uj, -q) + 

+C~"^/^g^ai(u;, q)ai{-uj, -q) - C^''^u;^ai(tJ, q)ai{-uj, -q) + C'^/'^a[{uj, q)a[{-uj, -q)+ (40) 
+2ibz*{qai{-uj, -q)a'Q{uj, q) +ujai{-uj, -9)02(0;, q) + E{uj, q)a[{-uj, -q))) , 



where we have omitted tildes for simplicity of notation and introduced the gauge-invariant electric 
field strength [33], 

E{uj, q) = uja2{u), q) + qao{uj, q) . (41) 
In addition we have Gauss's law ^ 



Together with 



coC^'/^a'oiu, q)+qC^l^a'^{vj, q) = Q 



E'{lo, q) = uja2{uj, q) + qa^iu}, q) 



eq. (42) gives 



;E', 



(42) 



(43) 



(44) 



This is an equation of motion for . To derive it replace 

0-2 a'2 — d2a^ , aj) — > tto — doa^ 
in tlie Lagrangian (|39|) and leave only terms linear in derivatives of , because only these will survive when we 



consider the equation of motion for in the = gauge. Then use the Fourier transform ( 37 1 
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4(0;, q) 



Plugging these expressions into the action ( 40 ) , we obtain 



S 



dz 



1 + h^z^ V C^/^ 



+ C,y^E^ + C,''l\'l + 2ibz\Eai)'^ . 



Corresponding fluctuation equations are 
2 / 1 



p. I - f - I o 1 1- (9/^)^(1 + + 

z Vi + ((i + ^^)^^)~^ Vi + &^2^ 1 -(g/^^)2(i + (i + 62)^4)-! 



+ UJ 



1 + (1 + 62)^4 



(l + 6224)(l + (l + ^2)_^4)3/2 



a'l+2z^ 



1 + 62 



2^2 



1 + (1 + 62)^4 1 + 62^4 



ai+ 



(45) 



(46) 



(47) 



1 + (1 + 62)^4 J ^ ' (1 + ^2^4)(l + (1 + 52)^4)1/2 

(48) 



_Z. Vanishing magnetic field 



In this subsection we set the magnetic field to zero. Fluctuations of E and ai fields then 
decouple, and we can consider separately transverse and longitudinal responses, 



E" + 



1 



z \1 + z- 



+ 2 1 



4^-2 



l-(g/a;)2(l + z4) 
l-(g/w)2(l + z4)-i 



S' + (a;2-g2(i + ^4)-i)^^0^ (49^ 



«'/ + 



2z-^ 



1 + z 



Oi + w 



1 + z^ 



ai = 0. 



(50) 



Let us first study the longitudinal response. In the near-horizon z^ 1 region eq. (49) becomes: 

(51) 



E" + -E' + uj^E = 0, 



The general solution of (51) is a linear combination of e^^'^^/z. We choose the solution with the 
incoming near-horizon behavior, since it corresponds to retarded propagator in the dual field theory 



E = C- 



(52) 



The constant C is undetermined, because the fluctuation equation is linear. When uz <^ 1, we 
obtain 



E = C{ - + iuj] . 



(53) 
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Condition (52) together with the boundary condition E(0) = (imposed to get normahzable 



solutions) defines an eigenvalue problem for the fluctuation equation (49). In the limit ujz,qz <C 1, 
(|49]) reduces to, 



E" + 



1 



z yi + z-"^ 

having as general solution, 



+ 2 1 



4\-2 



l-(g/a;)2(l + z4)-i 



E{z) = Ci + C2(g^ - 2u^)ViF (f sinh^i 



iz 



E' = 0, 



where F{z) is an elliptic integral of the first kind. In the limit z — t- oo it has an expansion 



^/iF [ i sinh 



-1 



iz 



1 



-K{l/2)+^+0^ ^ 



(54) 



(55) 



(56) 



where K{z) is the complete elliptic integral of the first kind. The solution (55) becomes in this 
limit 



E{z) = Ci- C2K{l/2){q^ - 



2C; 



2 2 
-ijj . 



Now we compare (53) and (57), and obtain as result 

((?2 - 2^2)^(1/2) 



lUJ 



2^2 



C. 



C2 



c 

'2^ 



(57) 



(58) 



Recalling the boundary condition E{{)) = 0, we deduce from (55) that Ci = 0, and consequently 

-1 



1 + 



ltd 



2uj^ 



K{l/2), 

which in the limit of small q, uj is solved by the considering leading orders in momentum g. 



iq- 



(59) 



(60) 



"^/2 4i^(l/2)- 

This excitation has been identified before, and is called [TGJ holographic zero sound. In the d = 2+1 
dimensional system this mode has been observed in [23J. Note that the speed of sound does not 
depend on dimensionality p of a probe brane and for any value of p is equal to the speed of the 
usual sound in the hydrodynamic regime [H]. In Section IV we will study current-current two- 
point functions, and the peak in the spectral function, corresponding to zero sound mode, will also 
be observed in the numerics. 



Now, let us consider the fluctuation equation (50) for the transverse gauge field component, in 



the limit cj, g <C 1. Then eq. (50) becomes 



2z^ 

ai + — — jOi = 0, 



1 + ^4 



(61) 
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with an exact solution being 

ai(z) = Ci + C2ViF (^ismh-\Viz)\ - l) . (62) 
In the near-horizon z — )• oo limit it is expanded as 

ai(z) ~Ci -C72K(l/2) + C2/z. (63) 



Comparing it with the incoming- wave solution (114), one obtains 



Ci = iiu; + K{l/2))C, C2 = C. (64) 



Then, near-boundary z <^ 1 expansion of ( 62 ) is given by 

ai{z) ^ A + Bz, (65) 

where 

A = {iuj + K{l/2))C, B = -C. (66) 

Therefore one may find the current two-point function {J^J^) = ^. In particular, its imaginary 
part is given by 



We provide numerical results for the transverse fluctuations in Section IV. 



2. Non-vanishing magnetic field 
In this subsection we are going to study the case of small magnetic field, 6^1, which will allow 



us to achieve some simplifications. Let us rewrite the action (46) as 



S = j dz{gEE''^ +UeE'^ + Qaa^ +Uaal+C^^\Eai)'^ , (68) 



where we have denoted 



. (l + (l + feV)^/^ (1 + (1 + 6V)V^ 

oT^w;:^ — — — ' ^ ^ 

^ (1 + (1 + 6V)V^ + a + 
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(1) 2ibz^ 

^ 1 + 62^4- ^'^) 



In the case of 6 <C 1, we can approximate 

,2 _ _3l_\ ' 1 + b'^z 



Qe = 7 ^ , = .,,24 ' C^^) 



(^ + -^)^^^ (^ + ^^^^^H-^-Tfe. ,,,, 

In the near-horizon hmit, for > 0, integrating the C^^^ term by parts, we arrive at 

Moreover, for z S> 1/Vb and z'^ S> l/(6((:<j2 — b'^q^Y^'^)-, we actually obtain decoupled system of 
equations 

E" - -E' + u'^E = ^ = ^1(1 - ia;z)e*'"^ , (76) 
z 



«i - -«i + (w^ - = ^ ai = ^2(1 - i V^2 _ 52^2^)gV'^'-'''9'^ . (77) 



z 

,2 ^ a2^2 



Now, assume that w ^ & y , and perform the linear transformation in (75) 



^ = «w(xi-^ + X2ai), ax = x\^-X2a\, (78) 
with arbitrary constant coefhcients xij X2) which brings the action to the form 

(79) 

Corresponding equations of motion are 

"^ + irT6^"^+" V' + ^(l + 62.4)J«i = 0- (81) 

The solutions are 

E = - ^-{\^iujz)e^^\ (82) 

Z UJ 
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ai 



z UJ 



-(1 =F iujz)e 



We impose the incoming-wave behavior, 

TP ■ ,'Xl+X2 



(Xi -X2)-(l-iwz) ) e" 



ai = [ ^^^^ + ixi + X2)-(l - lu^z) ] e-^ , 

Z Ul 



(83) 

(84) 
(85) 



which leaves us with two constant of integration xi =t X2- 



When cj ~ g ^ 1, we can consider fluctuation equations (49), (50), as for the case of vanishing 



magnetic field. Then we perform computations along the lines of the previous subsection, using 



now near- horizon boundary conditions (84) and (85) 



First, we match (84) in ujz <^ 1 limit, 



E = i(b{xi - X2) + iw^(x;i + X2)) + 



Xl + X2 



with eq. (57). Requiring that Ci = 0, we arrive at 



2 o 2 ^ - X2 I . 3 

q — 2(jj — — ujo h loj 

K{l/2)\ XI+X2 



0. 



Then, we match (85) in ujz <^ 1 limit. 



■ ( \ , f , ^b XI-X2 

01 = ^^{Xi - X2) + (Xi + X2)- H 

OJ z 



with eq. (63). Again, imposing normalizability condition Ci = 0, we obtain 

, , Xl + X2 , . 2 1 n 

w + I b \-iuj 1=0. 



K{l/2) \x1-X2 



Solving ( 87 ) together with ( 89 ) , we get 



262 



+ 



-(g2-4a;2) = 0. 



(86) 



(87) 



(88) 



(89) 



(90) 



[i^(l/2)]2 K{l/2) 

We see that in the presence of a magnetic field b zero sound mode develops a gap Uc in the spectrum, 

b 



K{l/2) ■ 



(91) 



Equivalently, we can obtain this result requiring that (87 1 and (891 have a non-trivial solution for xi iX2- 
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B. Effective theory for the sound mode 

Zero sound may also be studied in the framework of Ref. [37]. First, one introduces a hyper- 
surface z = 2;a in the bulk, integrating out degrees of freedom in the UV region < z < z\. The 
UV physics is then effectively encoded in the action by, 

S = ll d'x{fi{do(l) -Wo + wof - fi{d2ct> -W2 + W2f) , (92) 

where = 0^(2; = 0), Wfj, = a^{z = z\), and the "Godstone boson" (j) corresponds to breaking of 
the ?7(1) symmetry with a gauge field — Wfj_. The zero sound mode may be interpreted in such 
a framework as a mode coming from an excitation of the field (p, and therefore the speed of zero 
sound is given by the expression v = /2//o- Let us now compare the effective field theory action 
for the UV degrees of freedom with the bulk DBI action. To render the relation between bulk and 
boundary to be precise, we specify the zero boundary condition VF^ = 0, putting the Goldstone 
boson (j) to zero: 

S=ll Sx{flwl - flwl) (93) 
Let us consider all fields to be only z-dependent, in which case transverse fiuctuations decouple. 



and we can put these to zero. Then we can rewrite the bulk theory action (39) in a form 



where we have defined h{z) = y^l + (1 + b'^)z^. The solutions of the equations of motion on oq, 02, 
satisfying zero boundary condition at the AdS boundary, while being defined on the hypersurface 
z = ZA, are now given by: 

wo = Co / , W2 = C2 . (95) 







h^{z) ' Jq h{z 



To match the bulk action and the boundary theory ( 93 ) , we evaluate the action ( 94 ) on the solution 



of the EOM, which leaves us with the boundary terms aX z = z\ only 



S^\j d^'xiCowo - C2W2) , (96) 
which in turn with the help of ( |95| ), may be rewritten as ( |93| ) with 

2_ pdz(i+6V) 2_ r^ rf^(i+^v) 

"io h^z) ' "io h{z) ■ 

Therefore the speed of zero sound is given by 

2 _ dz{l + b'^z^) ( f'^ dz{l + b^z^)^ 
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When 6^1, one obtains 



1 

^ ^ 3r[l/4]2" 



■2- ^ , (99) 



and therefore for b'^z'^ <^ 1 one recovers the value of the speed of zero sound in vanishing b-field, 
uo = l/\/2, while for b'^z\ ^ 1 the speed of zero sound approaches zero. In this regime the 



description of the low energy physics by the effective action (92) presumably breaks down; it would 
be interesting to write the low energy description that would account for the gap in the spectrum. 



C. Thermodynamic properties of trivial embeddings 

We will study the thermodynamics of the trivial Dp brane embedding, to obtain as a result the 
value of the speed of the usual first (hydrodynamic) sound. We consider here the D3/Dp system 
with a 2 + 1 dimensional intersection, and in the Appendix we will study the supersymmetric 
D3/D7 system with a 3 + 1 dimensional intersection, in the presence of a non- vanishing magnetic 
field. 

The total prefactor of the action is irrelevant for the computation of the speed of first sound. 
The grand canonical potential is given by the equation 

E = -S= [ dp{p^ + &){p^ + B^ + d^r^'^ = a , (100) 



where a = T ( 1/4) /( 12-^7?). Using (35) one may calculate the charge density as. 



P = -^, (101) 

to find the energy density, being at zero temperature equal to the free energy, 

e = H + Ji^hP = MB"^ + d^^'^ . (102) 

Consequently, the speed of sound is given by 

P 9P ll + 2fe2 

de de 2 1 + ^ ^ 

Notice that this result is independent of p, which agrees with [35] . Observe that when the magnetic 

field vanishes we retrieve the value = 1/2, which we observed before in the dispersion relation 



(60). 



Notice also that all the steps performed in the above may be combined into one expression (use 

ds/dfich = d^y. 

2 dS/dfich Id log flch 



Pchd^s/dpl^ 2 aiogd ■ 



(104) 
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V. HOLOGRAPHIC CURRENT-CURRENT CORRELATORS AT FINITE 
FREQUENCY AND MOMENTUM 

In the previous section we have shown that a propagating mode (zero sound) develops a gap in 
the presence of the magnetic field. In this section we compute numerically the two-point function 
of the U{1) currents. First we set magnetic field to zero. We identify the holographic zero sound 
as a peak in the spectral function. We start by computing the density-density correlator (J^J^) 
using the linearized DBI action. We then proceed to computing the transverse correlator {J^J^). 
After that we proceed to the case of non-vanishing magnetic field and show that the gap in the 
zero sound spectrum shows itself on the numeric graphs. 

A. Fluctuations of electric field strength E 



Consider the fluctuation equation (47), near the boundary z = for any value of magnetic field: 

E" - {q^ - uj'^)E = 0. (105) 

Its general solution is of the form, 

E = AeFi + BeFii, (106) 
where we have denoted the two independent solutions as 

f, = l + 2i^.H(^!^.' + ..., (107) 

2 2 

Fii = z+ '^ ~ + • • • . (108) 
o 

The on-shell action is therefore given by 

du;dqAE{i^,q)AEi-u}, -q)^ n ' A z=e ■ (109) 

- u)^ AE{oJ,q) 
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Non-vanishing Green functions are 



{J'iu,q)j\-u;,-q)) = lim 



on— shell 



£-s>o 5ao{z = £, uj, q)5ao{z = e, —u, —q) 



on— shell 



6E{uj,q,z) 5E{—uj,—q,z 



lim 

6E{z = e, io, q)6E{z = e, —to, —q) dao{uj, q, z) 6ao{—uj, —q, z 
q^ BE{uJ,q) 
q^ - a;2 AEiyJ-.q) ' 



(J2(a;,g)j2(-a;,-(?)) = lim 



52 s-, 



on— shell 



e-s>0 ba2{^z = £, UJ, q)5a2{z = e, —oj, —q) 



on— shell 



6E{uj,q,z) 5E{—uj,—q,z 



lim 

e^o 6E{z = e, w, q)6E{z = e, -co, —q) 6a2{uj, q, z) 6a2{—uj, —q, z 
u? BE{uj,q) 
q^ - AE{uj,q) ' 



(J°(a;,g)j2(-a;,-(?)) = lim 



on— shell 



6a2{z = e, w, q)6ao{z = e, —uj, 

^'^ Son— shell 



6E{uj,q,z) 5E{—uj,—q,z 



lim 

e^o 6E{z = e, UJ, q)6E{z = e, —uj, —q) 6a2{uj, q, z) 6ao{—uj, —q, z 
ujq BE{uJ,q) 



(110) 



(111) 



(112) 



g2 - ^Eil^^q) ' 

Note that these expression agree with the Ward identity for the U{1) conserved current J^, 



uj{,f{uj, q)j\-uj, -q)) - q{j\io, q)j\-io, -q)) = , 



(113) 



We evaluate numerically the ratio Be/-Ae on the solution of equation (47) with incoming- wave 
near horizon behavior (52). In Fig. [3] we present numerical results for the real and imaginary 
parts of the Be/Ae for different values of w, g in the case of 6 = 0. The holographic zero sound 
corresponds to the peak in the spectral density. 



B. Fluctuations of the transverse component of the gauge field 

In this subsection we will compute numerically the holographic two-point function for the trans- 
verse current {J^{x)J^{y)). Let us put 6 = 0. 

In the near horizon regime z — )■ 00 the bulk solution, corresponding to the retarded current- 
current propagator in the dual field theory, takes the incoming-wave form 

ai = C— , (114) 

and in the vicinity of the boundary, the equation of motion becomes 

a'l - {q^ - a;2)ai = , (115) 
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FIG. 3: Real and imaginary parts of the Be/ Ae in the D3/Dp system with d = 2+1 dimensional intersection. 
The spectrum of excitations is exhausted by the holographic zero sound mode with the speed of sound 
uq ~ -j=, and the attenuation ~ g^. 



with a general solution being a combination of Fj and Fjj (107), (108) 



ai = AaFi + BaFii 



(116) 



Re[B/A] \ 




FIG. 4: Real and imaginary parts of the correlation function {J^{—q)J^{q)) in the D3/D7 system with d = 
2 + 1 dimensional intersection. No non-trivial collective excitation modes are observed. For small frequencies 
and momenta g ^ 1, the imaginary part of the correlation function behaves as lm[{J^{—q)J^{q))] ^ uj, 
independently of a particular value of q. 

The results of numerical evaluations of the holographic two-point function {J^{q)J^{—q)) = ^ 
are presented in Fig. [4| We see that it does not reveal any structure. 
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C. Non-vanishing magnetic field 

In the case of b ^ fluctuations of the longitudinal E{x^, x'^,z) and transverse ai{x^, x'^,z) 
components of the gauge field are no longer decoupled^. They are described by the action (68), 
which can be written as 

S = j dz (^(^-{QeE')' +UeE - ^(C(^))'ai^ E+ 

+ [GeEE' + Gaaia', + C(i)^ai]^=g° . (117) 

The first two lines vanish on shell. In the last line the cross term does not contribute to the 
variation of the on-shell action by the boundary z = values of the fields E and ai, because 

C«U=o = 0. (118) 

The on-shell action is then given by the boundary term 

Son-shell ~ hm / divdq ( — kEE' + aia[] . (119) 

Near the boundary the solutions to equations of motion are given by 

E = AeFi + BeFii , ai = AaFi + BaFu , (120) 



where Fj jj are defined by (107), (108) 



To compute current-current two-point function numerically, we follow [IH] , where general system 
of coupled equations in the bulk is studied. For arbitrary two independent solutions ^(i), ^(2) 



of the coupled system of fluctuation equations (47), (48), we define the matrix H = (^*(i), *^'(2))- 
Near the boundary it is expanded as 

H = AFi + BFii . (121) 



On-shell action (119) may be rewritten as 

J dujdq^'^ M^' , (122) 



^on— shell 



where 



1 



M = [ 9'^'^' ) . (123) 



We thank R. Davison for pointing this out to us. 
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The matrix of correlation functions is then given by (see eq. (2.34) in [46]) 



G ~ MBA-^ . 



(124) 



In such a form the current-current correlation matrix G is explicitly independent of a linear change 
of fields 



H ^ HR. 



(125) 



where R 



ri rs 
r2 r4 



is some arbitrary non-degenerate matrix. If $(1)^(2) 



£;(!), (2)' 

(i),{2) ) are some 
"1 



arbitrary independent solutions, then due to (|121|) we get 

A 



B 



(126) 



and therefore using (124) we obtain 

a 



V 



/I (2) _ aW 4(1)»(2) _ .(2) 



(127) 



Near-horizon solutions are given by (76), (77), which we can write as a linear combination of 
two independent solutions 



(1 - iujz)e'' 



(2) 



(1 - iujz)e'''^ 
-(1 - i^^uj"^ - &2g22)eV'^'-'''^ 



(128) 



Arbitrary near- horizon behavior, with the most general form (up to simultaneous rescaling of all 
fields by the same factor) may therefore be written as a linear combination of these two solutions, 



(1 - iujz)e^'^'' 



^^(1) + ^^(2) 



(129) 



On the other hand, fluctuation equations may be rewritten as 



J^i$" + f^2^'' + ^^$ = 0, 



(130) 



with matrices ill, 2, 3, being determined from (47), (48). Therefore linear combination of near 



horizon solutions (128) results in the same linear combination of the solutions near the boundary. 



Recall that the matrix correlation function (127) is the same for any such a non-degenerate linear 
combination. 



We therefore fix two arbitrary near- horizon conditions, say (128), determine corresponding 



coefficients A 



{1),(2) ^{1),{2) 



A^a'''''^' and B^a'^''^'^^ by integrating numerically fluctuation equations 
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FIG. 5: Real and imaginary parts of the component G22 of the correlation matrix ( 124) in the D3/D7 system 



with d — 2 + 1 dimensional intersection, for the magnetic field b — 0.001. The spectrum of excitations is 
exhausted by a gapped zero sound mode, with the value of the gap uic ~ b. 



(47), (48) up to the boundary and matching corresponding solutions with (120), and compute the 



correlation matrix (127). Each of the four components of the correlation matrix shows a gapped 
zero sound mode. 

In figure [5] we plot the real and imaginary parts of the G22 component, for b = 0.001. We see 
the gapped zero sound mode, with the gap which scales as Uc ~ b, in agreement with analytic 



result ujr 



K{l/2) 



of the previous section. 



VI. DISCUSSION 



In this paper we have studied current-current two-point functions at strong coupling. We have 
considered the current-current correlators at finite momenta, but did not observe any nontrivial 
structure in the spectral function, other than the zero sound^. 

It is instructive to compare the holographic density-density correlator with the form expected 
from the random phase approximation and reviewed in Section II. Within RPA the zero sound 
mode presents itself as a smeared delta- function like peak in Fig. 1, the Lindhard particle- hole 
continuum starts at q ~ w/vp and sharply ends at q ~ 2qi;'. The absence of the Lindhard 
continuum in the holographic computations can be explained by parametrically large values of the 



Landau parameters. The key point is eq. (11) which implies that since the zero sound velocity 



Note that our models are different from those studied in [71[S], where poles at finite momenta were observed in 
the holographic two-point functions of operators with nonvanishing charge under global (7(1). 
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that we observe is 0{1), the value of Fermi velocity scales like vp ^ l/\/FoFi. The regime of 
validity of our calculations is limited to w ~ q, and therefore the Lindhard continuum cannot be 
observed for parametrically large values of the Landau parameters. In the following we offer some 
speculations on how such a scenario can play out. 

We can argue that the Fermi velocity is parametrically small. Recall that q ~ q-v/A/^- Hence, 



eq. (60) implies that the zero sound attenuation is a ~ w'^^/X/fi. According to [47 \ this can be 



expressed in terms of the quasi-particle lifetime as 

a ~ - —VpFi ~ —FoFi (131) 
To derive the second approximate equality we used the Fermi liquid estimate 1/r ~ w^m*Fo^/q|. ~ 



w'^Fq / E*p. Eq. (131) implies that the Landau parameters are indeed parametrically large, -Fq-^I ~ 

We have analyzed the system in the presence of magnetic field and observed a gap in the 
excitation spectrum Wc. We derived a scaling relation Wc ~ B / ^. Note that the gap in the 
spectrum of non relativistic fermions scales linearly with while the relativistic fermions obey 
^/B scaling; Kohn's theorem implies that the gap in the spectrum of excitations is not changed 
when the pairwise interaction is turned on. In our setup charged fermions interact and can exchange 
momentum with A/" = 4 SYM degrees of freedom; the linear scaling of the gap with the magnetic 
field is consistent with the assumption that the effective degrees of freedom have an effective mass 



m* ~ ^. [According to eq. (10) this implies that Fi = 0(1); a scenario consistent with the 
discussion above may involve a parametrically large Fq ~ A^/^ but finite Fn, n > 0.] We do not 
quite understand the mechanism of dynamical mass generation at finite density ~ it is clearly very 
different from dynamical mass generation in a strongly interacting fermion system at zero density 

We already emphasized that a priori the very existence of zero sound is nontrivial, given the 
interaction of the charged matter with the uncharged superconformal degrees of freedom. It would 
be interesting to make this picture more precise and to see whether there is any relation to the 
recent studies of fermions in magnetic fields in the context of holography |48ft52j . It would also be 
interesting to compare our results with the correlators computed in the charged magnetic brane 
background [53] . 

At this point it is worth recalling the relation between the charge density and the value of the 



chemical potential, given by (35). As usual, the value of the charge density is proportional to d, 



The holographic dual of the latter involves repulsion of the probe brane from the bulk of of the AdS space; see e.g. 
1411 for a recent discussion. 
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p ~ Nc\^^~^^/'^d and the proportionality coefficient strongly depends on tfie dimensionality of the 
probe brane. The incompressibility dd/dfi is a smooth non-vanishing function of /i, b. This implies 
that we cannot rule out the existence of gapless modes in our system^^. Indeed, the analysis that 
led to the existence of the zero sound implicitly assumed a; ~ g ~ 6, and can be shown to break 
down for \uj\ < bq. We leave the search for gapless quasi normal modes for future work. Let us 
also note that a smooth compressibility is not compatible with the existence of Landau levels for 
the effective fermions. 

Appendix is devoted to the subject of higher derivative corrections to the DBI action for the 
probe Dp brane. The possibility of breakdown of the DBI description in the extreme infrared (very 
close to the horizon) was pointed out in [54J. Two possible causes were identified in the presence 
of the electric flux on the brane: strong back-reaction and vanishing of the effective string tension. 
The strength of back-reaction from the flavor branes is governed by the ratio Nf/Nc- We did not 
investigate 1/Nc corrections in the paper, although it is a very interesting problem. Instead, we 
explored the effects of the breakdown of the DBI description due to the vanishing of the effective 
string tension near the horizon. The strength of this effect is controlled by an inverse power of 
't Hooft coupling. In principle, such effects can be described by going to higher orders in the a' 
expansion of the effective action for open strings, which corresponds to adding higher derivative 
terms to the DBI Lagrangian. Unfortunately we are not aware of the precise structure of higher 
derivative corrections to DBI in the presence of the worldvolume electric field. However we were 
able to model this situation by writing generic higher derivative terms which become important 
near the horizon and completely change the effective metric for fluctuations there. 

The effect of such terms is confined to a very small region (which scales as an inverse power of 
't Hooft coupling in suitable units); outside of this region the second order differential equations 
derived from the DBI are applicable. In principle, one can solve the higher order fluctuation 
equation outwards from the horizon, and then feed the resulting solution into the second order 
equation. From the point of view of the latter, this amounts to modifying the boundary conditions: 
an outgoing wave (with a small coefficient) is added to the incoming wave near the horizon. We 
verify that this does not introduce any qualitative new features in the two-point functions. 



^ We thank D. Son for pointing this out to us. 
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Appendix: Higher-derivative corrections to L]jbi{o,i) 



The DBI description might break down in the near- horizon region |54j . and therefore higher 
derivative corrections become essential in that region. Consider higher derivative correction to the 
DBI Lagrangian of the form \55\ [56] 



e 



2v 55^^rt"''(V„FAp)(V;3F^.), (132) 
Alternatively, using the Bianchi identity, one may rewrite it as lyf^g^yy \F^^y aF"" . Here 

In this section we put L = 1. The induced AdS4, x metric on the trivially embedded Dp 
brane world- volume then takes the form 



ds" = p\-{dx^f + {dx^f + (dx^) + 4 + dnl, (133) 



This corresponds to non- vanishing Christoffel symbols in the AdS subspace, 



where ?7oo = —1) Vn = ^22 = 1- We fix the background value of Aq{p) (34) (with B = 0) and study 
the dynamics of the fiuctuation field ai(/0, x^^x"^). Consequently, the non- vanishing components of 
the field strength tensor covariant derivatives 

are given by 

ViFop = - Jfoi , VpFop = dpFop , (136) 
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ViFoi = p^^Fop , V2F01 = 52^01, VoFoi = 5oFoi - p^Fpi , VpFoi = ^pFoi - -Fqi , (137) 



V2-F12 — d2Fi2 + P^Fip , V0-F12 — 5oFi2 , VpFi2 — 9pi^l2 — , 



2 

p 



V2Fpi = SaFpi 



1 



21 



VoFpi — 9o-^pi -^01 ) VpFpi — 5pFpi . 



(138) 



(139) 



Let us now substitute the quantities (136)-(139) into the Lagrangian (132), which becomes in 
momentum representation, 



AL = e[-p\dpAoY-p\d'AoY + 



1 



2\2 



P p^\ p2 

4.o2„ ^2 , 4(a;2-g2) 



(140) 



To obtain the corrected equation of motion of the background field dpAQ, we put the ai fluctuations 
to zero and write the total, DBI + corrections, Lagrangian as 



L = p\l- (dpAoY - 4p\dpAoy + p'id'AoY 



The corresponding equation of motion 

p'dpAo 



is solved to first order in e by 

dpAo 



+ 2e[p^dpAo-dp{p''d^pAo)] = d 



n 



+ 5 do A 



(141) 



(142) 



(143) 



ip^ + S 

where we have denoted the correction to the background as 



5d.Ao 



(p4 + d^Y 



(144) 



Note that (144) approaches zero as 0{p^), near the horizon p = 0. Therefore the correction to the 
behavior of the background potential dpAQ does not substantially affect the near- horizon physics. 
Using the z = 1/p radial coordinate, and considering the near horizon limit ^ 1, we obtain 



from ( 140 ) the correction to the near-horizon DBI Lagrangian 

AL = e ((g2 - u;'')z\2af + (g^ - u;^)al) + {2a[ + za'lY) • 
This is to be added to the quadratic DBI near-horizon Lagrangian, 



(145) 



r 2/ /2 2 2\ 



(146) 
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As a result we obtain the following near-horizon Lagrangian: 

L = ((1 + 2e(g2 - a;2)) + 2e) + {-u? + e{q^ - u.'^f) z^a\ + 2e(za'2)' + ez2(a'/)' • (147) 

Up to a total derivative term^^ and 0{t) modification of the DBI behavior, this Lagrangian therefore 
may be rewritten as 

L = z\a'^ - uj\l) + ez\a'if , (148) 

with associated equation of motion 

a'i + ^a'l + w^a - e (^af + = . (149) 

To estimate the relative significance of the correction and DBI terms, let us compare terms a" and 
e (a'l" + fo'i"), when evaluated on the non-perturbed near-horizon solution e^^^ /z: 

"1 - ' e (^oi + -«i )-^—- (150) 

We observe that this correction is negligible. 

Unfortunately we are not aware of the exact form of the higher derivative corrections to the 
DBI action in the presence of the electric field on the world- volume of the probe Dp brane. In the 
following we will simply assume a particular expression for the higher derivative corrections to the 
Lagrangian for the transverse fluctuations: 

L = z^{dl - u?a\) + ez2+-(o;')2 ^ (151) 

with > 0. To estimate the significance of the correction term we need to compare contributions 
to the equation of motion from the terms 0(1) 

< = (152) 

and 0(e) 

e (z'^af + 2(1/ + 2)z'^-ia;" + (z^ + l)(i/ + 2)z^-'^a![) ~ ez^'-'^X + {ujzf) . (153) 

Therefore, if < u < 2, the correction becomes significant when z ^ (^tuflYiv (^^^ hierarchy 
of scales in Fig. [6|. If > 2, considering modes with sufficiently low frequency uj < e^/('^^2) the 
correction becomes significant when z 3> e^^^^"*^^ (see Fig. ItI). Finally, if > 2 and ui > €^^^'^~'^\ 



Corresponding boundary terms 2eza'i, evaluated on non-perturbed solution e"^^/z, vanish when z ^ 1. 
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— I 



FIG. 6: Hierarchy of scales in the near-horizon region for < < 2. 



FIG. 7: Hierarchy of scales in the near-horizon region for v > 2 and w < e 



l/{u-2) 



ion 



the Fig. jejis applicable, and the correction is significant when z » (^^^ly/v ■ Hence, in the regi 
z <C ^^J:^/^ the DBI description is vahd, provided that 0<i^<2ori/>2, uj > e^/^'^^'^\ The 
DBI description is valid in the region z <^ e^/i'^^'^) for > 2, u < e^/('^^^). 

The behavior of oi in the limit z ^ 1 where the DBI description is valid, is different from the 



incoming- wave (114): it has a qualitative form of "incoming wave" -|- 0{e) "outgoing wave". It 
is worth noting that the effect of higher derivative corrections on the current- current correlation 
function is essentially the same as an effect of non-zero 6-field. We verified that such a modification 
does not lead to any nontrivial structure in the spectral density. 
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